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Abstract 

In the framework of SU(2) gluo dynamics, we derive the tensor structure of the 
neutral gluon polarization tensor in an Abelian homogeneous magnetic field at finite 
temperature and calculate it in one-loop approximation in the Lorentz background 
field gauge. The imaginary time formalism and the Schwinger operator method are 
used. The latter is extended to the finite temperature case. The polarization tensor 
turns out to be non transversal. It can be written as a sum of ten tensor structures 
with corresponding form factors. Seven tensor structures are transversal, three are 
not. We represent the form factors in terms of double parametric integrals and the 
temperature sum which can be computed numerically. As applications we calculate 
the Debye mass and the magnetic mass of neutral gluons in the background field at 
high temperature. A comparison with the results of other authors is done. 

1 Introduction 



The investigations of QCD at high temperature carried out in recent years have elucidated 
the important role of color magnetic fields. In Refs. [1],[2] it was discovered in lattice 
simulations that sufficiently strong constant Abelian magnetic fields described by the 
potential of the form = BS^ 113 , where B is field strength, a is the index of internal 
symmetry, fi - Lorentz index, shift the deconfinement phase transition temperature T c . 
In particular, it was shown that an increase in the field strength decreases the transition 
temperature and for sufficiently strong field strengths T C (B) can be equal to zero. On 
the other hand, in Refs. [3] , [I] from the analysis of lattice simulations and in Refs. [6] , 
[7j, [8] from perturbative resummations of daisy graphs in the background field at high 
temperature it was found that Abelian chromo magnetic fields of order gB ~ g 4 T 2 , where 
g is a the gauge coupling, are spontaneously created. 
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These results are of interest not only for QCD but also for problems of the early 
universe where strong magnetic fields of different kind had likely been present [9]. They 
served as motivations for investigations began in our recent papers Refs.[10j. which 
goal is to determine the operator structure of the gluon polarization tensor in the con- 
stant Abelian chromomagnetic background field at finite temperature. This is necessary 
for investigations of the quark-gluon plasma (QGP), first of all, when resummations of 
perturbative series are carried out. As a preliminary step, the tensor structure of the 
gluon polarization tensor as well as the one-loop contributions to its form factors at zero 
temperature have been obtained and partially investigated therein. In the presence of the 
background field, for many reasons it is convenient to use the decomposition of the gauge 
fields in the internal space of the form = l/^/2(A 1 ± iA^^A^ = A^ and consider 
the former "charged" and the latter "neutral" gluons separately because of sufficiently 
different properties of them. That concerns not only physics but also the calculation 
procedures required for investigations. 

In the present paper in the framework of SU(2) gluodynamics, we derive the operator 
structure of the neutral gluon polarization tensor (PT) in the Abelian chromomagnetic 
background field at finite temperature. In actual calculations, as in Refs. |10| . [11] . 
we use the background Lorentz gauge and the Schwinger operator formalism based on 
the proper-time representation for propagators [12]. It is modified in order to account 
for the Matsubara imaginary time method at finite temperature and at the same time it 
preserves the Lorentz covariance when the momentum loop integrations are fulfilled. This 
generalizes the formulas derived earlier at zero temperature. As applications we calculate 
the one-loop contribution to the form factors. For instance, we derive the Debye mass and 
the "magnetic mass" of the neutral gluons in the background field at finite temperature. 
The limit of zero background field is also considered, to the correspondence with known 
results obtained already by other methods [15] is established. 

Here to present the results we would like to note that in the limit T — > oo the form 
factors can be calculated in terms of the Riemann Zeta-function. We also find that the 
"fictitious" pole (see for details, for instance, the surveys [15], [16]) in the gluon Green's 
function at finite temperature disappears in one-loop order if the external field is switched 
on. At the same time, the transversal neutral gluon field modes remain long range ones. 
They are not screened in the field at one-loop order, in contrast to the charged ones, as it 
was determined in Refs. [13], [8]. This difference is important. Its possible consequences 
will be discussed below in the main text. 

The paper is organized as follows. In the next section we introduce the necessary 
notation and review in brief the basic formulas. In sections 3, 4 we derive the operator 
structure of the gluon polarization tensor and develop the calculation procedure to carry 
out the integration over internal momenta of Feynman diagrams in the field at finite tem- 
perature. Explicit formulas for the form factors in the form of two-parametric integrals are 
obtained in section 5. In contrast to the previous paper [10], where the tadpole diagrams 
were not discussed in detail and instead some arguments in fewer of the cancellation of 
this contributions in the total by the surface terms appearing in other parts of the po- 
larization tensor were used. Here we consider this cancelation because of peculiarities 
appearing at finite temperature. In section 6 the transition to the zero temperature case 
is discussed. The Debye mass of the neutral gluon in the external field is calculated in 
section 7. In the next section we calculate the mean values of the operator in the physical 
states of the transverse modes and show that the gluon "magnetic mass" in the field is 
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zero in one-loop order although the fictitious pole of the Green function is eliminated. 
The discussion of the results obtained and further prospects are given in section 9. 

Throughout the paper we use latin letters a,b, . . . = 1, 2, 3 for the color indexes and 
Greek letters A,/i, ... = 1,...,4 for the Lorentz indices. Summation over doubly ap- 
pearing indices is assumed. All formulas are in the Euclidean formulation. We put all 
constants including the coupling equal to unity. Since the present work is a continuation 
of investigations began in Ref. [10J, we follow the notations, definitions and calculation 
procedures used therein as close as possible. 



2 Basic Notations 

In this section we collect the well known basic formulas for SU(2) gluodynamics to set up 
the notations which we will use. We work in the Euclidean version. Dropping arguments 
and indices, the Lagrangian is 

£ = ~F% ~ ^(d,A,f + r) d,D, V , (1) 

where £ is the gauge fixing parameter and rj is the ghost field. The action is S — J dx L 
and the generating functional of the Green functions is Z — J DAexp(S). In the following 
we divide the gauge field A*(x) into background field B°(x) and quantum fluctuations 

Oft*), 

A^(x)=B;(x)+Q;(x). (2) 
The covariant derivative depending on a field A is 

d_ 

and the field strength is 

JLa'(x) - — 
dx» uK } dx" 

and 



Da M = ^ ab + ^ acbA » ( 3 ) 



F ^ = TS^C*) - T^ A » + < abCA » A ^) ( 4 ) 



[D,[AlD v [A\\ ah = e acb F^[A] (5) 
holds. For the field splited into background and quantum parts we note 

F; u [B + Q] = F^[B] + Df[B]Q b v {x)-Dl b [B]Ql{x) 

+e abc Q b l (x)Ql(x). (6) 

The square of it is 

-\(F«AB + Q]) 2 = -\{f; v [B]) 2 + QlD^[B]Fl v [B] 

- l -Q a ,Kf u Q b v + M, + M,. (7) 

The second term in the r.h.s. is linear in the quantum field and disappears if the back- 
ground fulfills its equation of motion which will hold in our case of a constant background 
field. The third term is quadratic in and it defines the 'free part' with the kernel 

K* = -5, v D a ;[B]Df[B] + D^[B]D?[B] - 2e acb F^[B]. (8) 



3 



The interaction of the quantum field is represented by the vertex factors 

M 3 = -e abc (D; d Qt) CftQl, 

m 4 = -)qiqiq\qI + \qiqIqiqI. (9) 



The complete Lagrangian 



C = -l(F;jB + Q}) 2 + £ gi + £ gh (10) 



consists of ([7]), the gauge fixing term (in the following we put £ = 1), 

£gf = ~ {D;[B]Ql) 2 = ^Q a »D a ;[B]D*[B]Ql (11) 

and the ghost term 

C gh = rfD; c [B] {pf[B\ + e"°Qj) rj b . (12) 

These formulas are valid for an arbitrary background field. Now we turn to the specific 
background of an Abelian homogeneous magnetic field of strength B which is oriented 
along the third axis in both, color and configuration space. An explicit representation of 
its vector potential is 

and the corresponding field strength is 



B*Jx) = S a3 6^x 2 B (13) 



Ff. = S^Fij = Be 3ij , (14) 

where only the spatial components = 1,2,3) are nonzero. Once the background is 
chosen Abelian it is useful to turn into the so called charged basis, 

Qit = Ql (is) 

with the interpretation of as color charged fields and as color neutral field. This 
is in parallel to electrically charged and neutral fields. Note also that is real while 
are complex conjugated one to the other. In the following we will omit the word color 
when speaking about charged and neutral objects. The same transformation is done for 
the ghosts, 

V, = Vl (16) 
A summation over the color indices turns into 

QlQl = Q,Qu + w+w- + w~w+. (17) 

All appearing quantities have to be transformed into that basis. For the covariant deriva- 
tive we obtain 

Df = d„ D~+ = 8»-iB, = D, 

(18) 

D+- = d^ + iB, ee D* 



where D* is the complex conjugated to D^. Starting from here we do not need any longer 
to indicate the arguments in the covariant derivatives. 

Before proceeding with writing down the remaining formulas in the charged basis it 
is useful to turn into momentum representation. This can be done in a standard way by 
the formal rules. It remains to define the signs in the exponential factors. We adopt the 
notation 

Q ~ e~ lkx , W~ ~ e- ipx , W + ~ e ip ' x . (19) 

In all following calculation the momentum k will denote the momentum of a neutral line 
and the momenta p and p' that of the charged lines whereby k and p are incoming and p' 
is outgoing. In these notations the covariant derivative acts on a W~ and turns into 

D M = -itidp + = -ijy (20) 

Note that the components of the momentum p^ do not commute, 

\p ttlPv ] = iBF tlv . (21) 

In these notations the quadratic term of the action turns into 

-\QlKQl = (22) 

\QpK%Q v + \w+k;+w- + \w;k+~w: 



with 
and 



^ iiu — K^yifc) — 5^ v k k^k u (23) 



= Kp V (p) = 5^p 2 - Pp Pv + 2iBF fiv . (24) 

We use the arguments k and p instead of the indices to indicate to which line a K^ u 
belongs. The third term in the r.h.s. of Eq. (l22|) is the same as the second one due to the 
complex conjugation rules. In the Feynman rules (i^) -1 is the line for neutral gluons 
and is denoted by a wavy line and (i^"^) -1 is the line for charged gluons and is denoted by 
a directed solid line. We remark that these lines represent propagators in the background 
of the magnetic field. Frequently they are denoted by thick or double lines. Because we 
have in this paper no other lines the notation with ordinary (thin) lines is unique. 
Here we note that the spectrum of the operator fl24l in a constant magnetic field, 



E 2 n = p\ + B{2n + 1), n = -1, 0, 1, (25) 

contains a tachyonic mode at n — — 1. p% is momentum along the field direction B = B 3 . 
This state is a peculiar of non Abelian gauge fields as it is discussed in different aspects 
in the literature (see, for instance Ref.pJ]],|7| and references therein). 

For later use we introduce here the set of eigenstates for the operator (123]) . For 
the color neutral states we take exactly the same polarizations | k, s > as known from 
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Figure 1: Notations for the vertex .M3 
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(26) 



with h 



k\ ~\~ ^ k 



ftf 



b 2 -r fcf. Here the polarizations s — 1,2 describe the two 
transversal gluons (k^ \ k, s — 1, 2 > M = 0), s = 3 is the longitudinal one and s = 4 after 
rotation into Minkowski space becomes the time like one. For the transversal gluons 



K^(k) I k,s = 1,2 > u = (kl + k 2 ) |M = 1,2>, 



(27) 



holds. 

After discussing the free part of the Lagrangian (flQj) . (|7j), in the 'charged basis' we 
note that for the vertex factor Ai 3 in (Q we obtain 



M 3 = W-T^xW+Q 



:28) 



with 



r, 



5„ 



uv\ = <v( fc - 2p) x + 5 X n(p + k) v + 5 Xv (p - 2fc) M . (29) 
The notations are shown in Fig{TJ It should be remarked that all graphs and combinatorial 
factors are exactly the same as in the well known case without magnetic field. On this 
level the only difference is in the meaning of the momentum p^ which in our case depends 
on the background magnetic field, see Eq. (l20|) . 

The vertex factor describing the interaction of the neutral gluons with charged ghost 
fields is 

Mf = v-rfv+Qx (so) 

where 

TT=p'x = (p + k)x. (31) 

We also need in the four particle vertexes which are momentum independent and have 
the same form as at zero external field. 
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3 Operator Structures 



The neutral PT is denoted by U xx > (k) where the argument k is an ordinary momentum. Its 
one-loop diagram representation is shown in Fig. [2J In this section we discuss the general 
tensor structure of it at zero and finite temperature. As it was shown in Refs. [TU].|llj 
it is not transversal in a magnetic background field. That means that the condition 
k\H\\r(k) = does not hold. This follows either from the Slavnov- Taylor identity for the 
gluon Green function or from an explicit one-loop calculation. So, we are left with the 
weaker condition 

k x U X y(k)ky = 0. (32) 

At zero temperature, it can be combined with the remaining Lorentz symmetry which 
results in a dependence of U X y(k) on two vectors, l x and h\, and on the magnetic field. 
We use the notations 



( \ 
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\ h J 

( 
-l 

o 

V o 



f h\ 
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V o j 



( h \ 

-h 






10 0^ 





/ 



(33) 



The third vector is d, 



F k 



Note that here and further in actual calculations the 
magnetic field strength B is put equal to unity. For the vectors k x = l x + h x holds. 

The general structure of U X y(k) at T = allowed by (|32|) and the vectors l x and h x 
is determined by the set of tensor structures 



XX' 


— l 2 $ xx > — hh' 




T (2) 
XX' 


= h 2 5 xx , - h x h y 


= d x dy 


T (3) 
XX' 


= h^ly+l^y 


- hhy - 


AA' 


= i(l X dy - d X ly) 


+ il 2 F xx 


T (5) 
AA' 


= hH\ xl -Hi xl 




T (6) 
AA' 


= iF X y 





I'Xi'X' 



(34) 



together with the identity d x hy — h x d 



y /iA 
^ XX' 



diag(l, 1, 0, 0) and S 



T {i) K, 
1 XX' K x' 



fiX 

with i 



= h F X y. Further we introduced the notations 
diag(0,0, 1, 1). The first four operators are transversal, 

1,2,3,4, the last two fulfill (I3"2l . only. The sum of the 



first three operators is just the transversal part of the kernel of the quadratic part of the 
action, Eq. (|23|) . 



Ti\]+nt>+T : 



(2) 



,(3) 



AA 



AA' 



. Xx ,-K xx ,(k). (35) 

At finite temperature, an additional vector - the velocity of the thermostat - must be 
accounted for and used in construction of the tensors . Therefore new tensor structures 
appear. We chose them in the form: 



T 



(7) 



AA' 



{u X ly + l X Uy){uk) - 5 x \,(uk) 2 - U X Uyf 
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(36) 



Obviously the sum 




(u X ky + k X Uy)(uk) - 5 X \'(uk) 2 - U x Uyk 2 = B\\i 



(37) 



equals to one of two transversal tensor structures commonly used at zero field [15] . In 
that case the first structure is given by the sum of tensors Eqs. (|35|) . (1231) . Below in 
actual calculations we use the reference frame of the thermostat, so only one component 
of is nonzero: = (0, 0, 0, 1)^. We mention that the first two tensors in Eq. (l36|) are 
transversal and the other two satisfy the weaker condition (13"!?|) . 

We adopt the following way for the representation of our expressions. The dimen- 
sionality of the polarization tensor is implemented in the tensors Tw. To restore the 
dimensionality for the tensors in Eqs. (1341) and fl36l) . one has to multiply the operator 
by the factor B, and the operator by yB. The form factors are dimension less 
functions of dimension less momenta I 2 , h 2 , and temperature T. That means, in fact they 
are measured in units of B. To restore the correct dimensionality one has to replace 
I 2 — > l 2 /B, h 2 — > h 2 /B, and T — > T/\/B. Correspondingly, the arguments of all functions 
appearing in the actual calculations are also dimension less. 

Knowing the operators (I34j) and (136]) . which may appear at zero and finite temperature, 
the polarization tensor can be represented in the form 



where the form factors U^(k) depend on the external momentum k^ through the variables 
I 2 and h 2 at zero temperature and h 2 , k^ and k^ at finite temperature. In the former case, 
the polarization tensor U\y(k) is real and symmetric in its indices, so the form factors 
n( 4 )(fc) and n( 6 )(/c) are zero. In the latter case all structures will contribute in general. 

In the Matsubara formalism, it is possible to add to the set of tensors Eq. (1361) a 
structure of a special type which contributes in the static limit, k^ = 0, only (remember, 
fc 4 takes discrete values). It should be of the form u^u^iki) with ${k^) possesses the 
following property: it is non zero for fc 4 = 0, only. This term is obviously transversal for 
itself. So, the set of tensors in Eq. (l38|) could be extended due to this structure. On the 
other side this structure is a linear combination of T$, and T®, at k± = 0. 

4 Calculation of the Neutral Polarization Tensor 

In this section we calculate the one-loop contribution to the PT at finite temperature. 
The imaginary time formalism is used. That means that in loops the integration over 
momentum component p4 is substituted by an infinite series in discrete values = 2tcNT: 

+oo +oo 



n AA ,(A0 = E n«(A;)T ; 



(38) 



8=1 




/ $/(p 4 )-T E f(2*NT). 
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The neutral polarization tensor has the following representation in momentum space 
^see Fig. |2D 



U xy (k)=T J Il(p,p 4 ,k,k 4/ 



, (27r)3 "v— --^v, (39) 

where in the integrand we noted explicitly the dependence on the external momentum 
and the momentum inside loops. In what follows, as in Ref . [TU] , at intermediate stage of 
computation we, for brevity, shall omit the general factors and the signs of integration 
and summation. That is, we relate these factors and operations with the momentum p 
standing in r.h.s. of equations. Within this convention we write 

n A A'(/c) = ^ nv\G nii'{p)^ ii'v'\'Gv>v{p — k) (40) 
-P\G(j>)(p - k) X 'G(p -k) - (p- k) x G(p)px'G(p - k) 
+G xy (p) + G x , x {p) - 25 xx ,TrG(p). 

The second line gives the contribution from the ghost loops and the third one is due to 
the tadpole diagram. 





Figure 2: The neutral polarization tensor 



The vertex factor is given in Eq. (1291) . For a convenient grouping of terms it is useful 
to rearrange it, 

-p 

g ilv {k-2p) x + 2{g Xil k v - g Xv k^ + g XfJt (p - k) v + gxuPn, 



= r (i) , r (2) , r (3) 

where in the last line a subdivision into three parts is done. 
The propagators are given by 

1 f°° 2 

G(p) = -,= \ ds e~ sp 



p2 



„2 



(p-k) 

for the scalar lines and by 



G X y{p) 



(p — k) 2 + 2iF / 

[°°dt e-^-^EZ* 
Jo 



(41) 



G{p -k)= 1 = I" dt e~ t{ - p - k) " (42) 



GxAp-k) = [ - -i ■ — 1 (43) 



for the vector lines (in the Feynman gauge, £ = 1) with 

El, = (e 2 -) AA , 

= Six + iF w sinh(2s) + 8^ x , cosh(2s). (44) 

At zero temperature, the momentum integration can be carried out by means of 
Schwinger's algebraic procedure [12] and converted into an integration over two scalar 
parameters, s and t. Here we educe the known results in order to present their modifica- 
tions at T ^ 0. The basic exponential is 

6 = e -^ e -tiv-k? ( 45 ) 

and the integration over the momentum p is denoted by the average (•••)• The following 
formulas hold: 

exp \-k + &S X ) k 

(Q) = —± ^ S+T 1 J (46) 

V 1 (47r) 2 (s + t)sinh(s + i) 1 ' 

with S = tanh(s) and T = tanh(t) and 



(P^e> = (p) (e>, (47) 
(p^e> = ff4*l (4*) W (48) 



D J^\D ) v KD 1 / iaJ/ 



The notation A = E l — \ and D = E s+t — 1 is used. Explicit formulas are 



A t .^sinh(s) sinh(t) ^ ± cosh(s) sinh(t) 

D s + t sinh(s + 1) sinh(s + 1) 



(49) 



along with 



-1iFE- s 5u smh(s-t) P cosh(s-t) , , 

— = — iF — r) z + $± — 77 f, (50) 

D T s + t smh(s + t) smh(s + t) 

where we dropped the indices. It should be remarked that all these matrices, i.e., E, F, 
D and A commute. In addition we need the relation 

pis), = e- sp2 Pfl e sp2 = E\ vVv 

for commuting a factor p, with the propagator G(p), 

PtiG^'ip) =G(p)p fl >. (51) 

Now we turn to the finite temperature case. Our goal is to account for the temperature 
dependence within the above representation in a natural way. Usually in the imaginary 
time formalism the summation over p 4 and the integration over three-momenta are carried 
out separately. To restore the equivalence of these variables and to make use of the 
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formulas (|46p - (|48p we proceed in the following way. First we note that any function 
/(p4 = 2ft NT) of pi can be written in the form 



/(p 4 = 2nNT) = J dpJ( Pi )5( Pi - 2ixNT) 

+00 

= / dp 4 f{p,)— J dXe^~^ NT \ (52) 

— 00 

where S(x) is Dirac's delta-function. Then we change the order of integration in variables 
A and p± and fulfill the momentum integration as at zero temperature. The factor l/(2ir) 
coming from the delta-function and the factor l/(27r) 3 coming from the three-momentum 
integration give in the product the factor l/(27r) 4 appearing at zero temperature. The only 
new factor, e iXp4 , appears as another exponential in Eq. (145 jl . Since the fourth component 
is not related with the magnetic field, the calculation is actually the same as at zero 
temperature. Remind that in course of it the factor l/(27r) 4 results in the factor l/(47r) 2 
in the function (Q(s,t)) [12]. The integral over p 4 is Gaussian, and we obtain for the 
algebraic averaging (or bracket) procedure (...) at A 7^ 

(Q(s,t, X)) x = (e iA ^6M)> = (©(M)> exp L^ + iXkJ-j , (53) 

where we marked the averaging procedure with A-dependence by the subscript A. Here k± 
is the discrete fourth component of the external momentum and we introduced q = s + 1 
as a convenient variable. The expression in the angle brackets in the r.h.s. is the zero 
temperature value given in Eq. PSl) . Below we will denote the function f[53"j) as Q(s,t, A). 

At the next step, we integrate the function 0(s,t, A) over A and after the restoration 
of the sum over N we obtain the basic expression at finite temperature: 



+00 +9° 

(Q{s,t)) T = T £ / dXQ{s,t,X)e- iX(2nNT) 



+00 



= T E V 47r 9( ( s > *)> ex P ~ k l~ + 2ht(27tNT) - (2nNT) 2 q . (54) 

We denote the function standing under the sign by 0t- 

To obtain the result of the bracket procedure with the momentum p^ entering, one has 
to differentiate Eq.( l53l) with respect to iX and then to calculate the integral over A. This 
can be done also by means of differentiation of Eq.( l54l) with respect to the parameter 
6jv = 2k NT. In this way the all integrals of interest can be computed. 

The expressions with the spatial indexes i, j = 1, 2, 3 are given by Eqs. fl46l) - fj48l) . where 
the function (|54"|) must be substituted. 

The average with one spatial component and p± is 

<pa*e(8,t)) T = t g 7 dXe ^NT) k. + iX/2 ,a\ 

+°° /A \ 

= E \n k ). 2irNTQ(s,t,X). (55) 



N=-oo 
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The bracket procedure with p\ results in the following expression: 



9 / 2 <^ 

E tM) (27riVT) 2 . (56) 



N=-oo 



Note as an interesting fact that the /^-dependence in the last two equations comes in 
through the exponential, only. 

Thus, we collected the necessary integrals over internal momentum which appear in 
the magnetic background field at finite temperature. These expressions are useful when 
the high temperature limit is investigated because due to the exponential factor in Eq. (154p 
a few first terms in the series (N = 0,1,2,...) give the leading contributions at T — > oo. 
The factor yjq entering the integrals over s and t parameters ensures the convergence at 
s,t — > 0. This corresponds to the superficial divergence degree. In the four dimensional 
theory the form factors have zero superficial divergence degree, hence the three dimen- 
sional theory which effectively appears in the high temperature expansion is ultraviolet 
finite. The convergence at infinity is due to the exponential factors except for the com- 
plications resulting from the tachyonic mode which will be discussed later. However, the 
low temperature limit is less trivial. 

Now we consider another representation for the integrals, convenient at low tempera- 
ture and carry our the ultraviolet renormalization. For this purpose we make resumma- 
tions of the form 



N=-oo " Z TV 



E e^ N2+ ^ = J- E (57) 



of the series in N in the Eq. (154|) with the parameters a = AnTk^t and z = 4n 2 T 2 q 
chosen. In this case the dependence on k\ in the exponential disappears and we obtain 
for (Q(s,t)) T : 

+oc ( N 2 kJN\ 

<e(M)) r = E <©M)> exp -—+*-*— . (58) 



N=—oo 



AT 2 q qT 



This is just the representation of interest. As above, the function standing under the sign 
of the sum will be denoted by 0^- 

To obtain the functions (0(s, t)p±)T an d t)pl)r we have to calculate the integrals 
over the variable A with the powers of A and A 2 in the Eqs.f l55j) and (1561) by means of 
differentiation of the expression (I58p with respect to the parameter = ik^t/q one and 
two times, correspondingly. In this way we obtain in the first case, 

/ ^/ ^ ^ (, t iN \ ^ ( N2 ihNt\ 

ws)*-^ + < e (^» + (59) 

This expression can be combined with the bracket operation for the spatial momenta p^. 
In fact, the function (4^)^ for [i = 4 equals to (t/q)k 4: that coincides with the first term 
in the bracket in the above equation. To account for the second term, we assume that 
the thermostat is at rest and therefore the vector w M has only one non zero component, 
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v>n = 5^4- To write down Eq. fjoUj) in short, we introduce the notation = ytU^. Then, 
the average operation with one internal momentum yields 

+00 A 

(Pm©(*,*)>t= E + «M)Q(*.*)r. (60) 

N=— 00 

In the same manner, in case of two internal momenta the result can be presented in the 
form, 



N=-oo 



e T (s,t). (6i) 



Thus, we have calculated all integrals needed for what follows. 

The derived expressions ( H6l) -( )48l) and ( 1581) . ( 1601) . (I6T1) solve the problem on the momen- 
tum integration in the Schwinger operator formalism at finite temperature. The results 
are expressed in terms of two-parametric integrals, as at zero temperature, and a sum 
over N which takes into account the temperature dependence. 

So, all further steps of calculations necessary to obtain the expression for the polariza- 
tion tensor actually coincide with that in Ref.[10j. Within this formalism, the polarization 
tensor becomes an expression of the type 



n A A'(£0=/ / dsdt(M xy (p,k)e) Tj (62) 
Jo Jo 



where in M xx >(p,k) we collected all factors appearing from the vertexes and from the 
lines except for that which go into 0. By using fl60l) and fIBTl) the average (the momentum 
integration over p) can be transformed into 

(M xx ,(p,k)0) T = M xx ,(s,t)0 T , (63) 

where now M x y(s, t) collects all factors except for Q T . 

As in Ref.|10j. we break the whole polarization tensor into parts according to the 
division introduced in Eq. fflTJ) . We write 

Hay (k) = E n& (*) + nf A ? st (*) + Tl xx , (k)^ (64) 
hi 

with i,j = 1, 2, 3 and 

= ^xG^{p)Y% x G u , u {p - k) (65) 

including corresponding subdivisions of M. 

To organize our calculation, we remind that the PT is calculated in terms of double 
parametric integrals in s and t. The function (0(s, £)) fj46|) is symmetric whereas 0y(s, t) 

is not, because of the factor e ( i T ' = ^\ uk ) 1 in the temperature dependent exponential. 
To restore the s «-> t symmetry, we write 

e 2(uk)t = St + At ( 66 ) 

where 

S T = - (e^ + e 2 ^ s ) , A T = - (e 2 ^ - e 2 ^ s ) . (67) 
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Then the function Qx(s,t) can be split in symmetric and antisymmetric parts, ©7. = 
Of. + Of, with 



Of, = (G(s,t))S T e-^, 

_ N 2 

Gf, = (e(s,t))A T e (68) 

With these definitions introduced, the terms entering the integral in Eq. (1621) have the 
following general structure. The expressions in M\y(s, t) which are symmetric under s <-> 
go multiplied by the symmetric temperature factor Of,, and vise versa, the antisymmetric 
functions go multiplied by Of. This observation gives us the possibility to make use of 
the results in Ref.[TU], where the symmetric form factors for the operators — 
where calculated. At finite temperature, they should be multiplied by the factor Of. in 
the total expression. The antisymmetric terms as well as the terms giving contributions 
to the remaining form factors must be calculated separately. 

In fact, this procedure is rather simple and it is described in the previous work. To 
make the present paper self-contained, we educe the details of this calculation in the 
Appendix. In the main text, we restrict ourselves to the description of it in general 
words. For the first six operators, i.e., for that which do not contain the vectors u^, the 
terms appearing after the momentum integrations have the following forms. First, 

ITaa' = PxPy + aS l l y + bSiy + icF X y , (69) 

where Pa is given in terms of the vectors (1331) . 

P x = r l x + a id x + (3 h x (70) 

and Py is the transposed expression, r, a, /3, a, b, c are some functions specific for different 
parts of IIaa'- 

The second type of expressions has a slightly more complicated form 

IIaa' = PxQx> + QjPy + aS l l x , + b5i x , + icFxx, (71) 
with Pa from (1A.2I) and 

Qx = s / A + 7 idx + 5 h x , (72) 

where s, 7, 6 are some other functions. Then, from the requirement that the weak transver- 
sality condition fl32|) holds, one can derive the expressions standing in front of the operators 
TW — This procedure is efficient but, of course, not obligate. 

The form factors standing in front of the operators — can be obtained after 
integration in accordance with Eqs. fl60l) . fl6~T|) as that part of the expressions which are 
proportional to one or to two powers of u^. To make clear the structure of the expressions 
which appear after the averaging procedure, we note that it results in the formal substi- 
tution of in the initial expressions by p M — >• + = [(jjjk) + u]^ in the final ones. So, 
the M-dependent parts are easily determined. The necessary details for different parts of 
the polarization tensor are given in the next section. 

We complete this section with the description of the renormalization procedure adopted. 
As it is well known, the divergent parts of the polarization tensor do not depend on the 
temperature and field. So, each form factor can be written as follows 

n (i) (p,T) = [n (i) (p,T)-n (i) (p = o,T = o)] + n (i) (p = o,T = o). (73) 
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Then, the expression in the brackets is finite whereas the last term is divergent and 
must be renormalized by the standard procedure in quantum field theory. In terms of 
the resummed series ( l58|) . (!60!) . (!6l!) the term with N = just corresponds to the zero 
temperature case. So, actually the above procedure has relevance to these terms, only. 
As a result, we obtain the renormalized polarization tensor at finite temperature in the 
field presence which is the object of interest. 



5 Calculation of the form factors 



In this section we calculate the contributions to the form factors stemming from individual 
= 1,2,3) introduced in Eq. (l4Tj) . At zero temperature that was done in 



terms T 



Ref.[10j. So, here we present mainly that part of calculations having relevance to the 
temperature dependence of these expressions. The first is 



n 



ii 

AA' 



(k - 2p) x G^(p)(k - 2p) A ,GU,(p - k) 



(74) 



and it transforms into 



Mft(M)©r = <(* - 2^)a(A: - 2p(s))yE;°,E^ 6 



r 



'D 



k - 2u 



-Ai 



E~ S F 



trE' 



-s-t 



AA' 



1 - 



2E S 



A 
D 



k-2u 



(75) 



We note the property E s ^ = (4^) . The trace is 

trE' 8 - 1 = 2 [1 + cosh(2g)] . 



(76) 



Remind the variable q = s + t. 

Here we also introduce the variable £ = s — t which is antisymmetric with respect to 
the replacement s <->• t and assume s and t to be replaced, s = (g + £)/2 and t — (q — £)/2. 

Using the notation of Eq. flA.lj) we define 



Pa 



1 



2±U 



£ Bsinh(t) sinh(£) 

-l x + 2id x —rr + h x - 



q sinh(g) 

rl x + iad x + (3h x 



sinh(g) 



(77) 



and below we will express the obtained results for different form factors in terms of r, a, (3. 

The second part in the expression f!75|) has to be integrated by parts. This procedure 
at finite temperature requires a special explanation. We represent 



-2iFE~ 



'rS"- i F C -^ + S^ + C 
smh(g) 



(78) 



where C is a constant, i.e., it must be independent on £. We include it in the integration 
procedure in order to make use of this parameter, see below. The derivative with respect 
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to £ will be integrated by parts. We should note that expressions which are symmetric 
under an exchange of s and t, i.e., which depend on q only, are not affected in Eq. (!78l . 
So, in 0t we have to differentiate only the terms which depend on £, 

JLe T =C-B l - 2(uk)) T , (79) 



d£ \2 
with the notation 

B 1 = rl 2 + f3h 2 . (80) 
Then, after integration by parts, the last term in the Eq. (!75l) gives the contribution 



' E~ s F\ 

- 4i (-^f-J = 2(1 + cosh(2q))(2(uk) - Bx) 

smh(q) 



that has to be added to the first part of the equation. It contains a term Au^u^Qt which 
is is the quadratic in w M . Now we chose the constant C = —2j^ in a way that this term 
cancels and considerably simplifications in further calculations appear. In this a way we 
arrive at the final expression at finite temperature. As concerns the surface term, it is 
canceled by the contribution of the tadpole diagram, see below. 

Applying formula (IA.3I) from the Appendix for the function flA.ll) with the parameters 
a = rBi and b = (3Bi , giving the contributions to the first six form factors, and gathering 
the factors at the -^-dependent structures giving rise to the operators — T( 10 \ we obtain 
for Ml\,: 

M n ( S ,t) = {- r 2 TW + (a 2 -p 2 )TW-rf3T® 
-arT (4) + [B x coth(g) - 2{uk)}T {6} 

+2{uk)[-rT i7) - /5T (8) + aT (9) ]} 2(1 + cosh(2g)). (82) 

This rather simple expression includes the symmetric and antisymmetric with respect to 
the t terms. Since the integral in s, t is symmetric, actually, the factor 0^ stands at 
the former terms, and 0^ at the latter ones. This remark concerns all the expressions 
written below for other parts of H\\'. 

The next contribution is M 22 . From (HT|) we get 



n AV = 4 (S^kv - 5 X uk^) G flfl f(p) {5 X ' t i'k vl - <WV) G vv{v ~ k) (83) 



and 



Afg.(V) = i(E-^(kE- l k)+E{ y (kE-k) (84) 
-(E-'k) x (E-<k) y -(E'k) x (E-k h ,). 

Here we have expression of the type of (IA.4|) . The parameters are 



r = 1, a = sinh(2t), (3 = cosh(2t), 
s — 1, 7 = sinh(2s), 5 = cosh(2s), 



(85) 
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and 



a = -2l 2 - /i 2 (cosh(2s) + cosh(2t)), 

b = -Z 2 (cosh(2s) + cosh(2t)) - 2h 2 cosh(2s) cosh(2t). 



Using formula f]A.6|) from the Appendix we obtain 



M 22 {s,t) = 8T« + 8cosh(2(s + t))T( 2 ) 
+4(cosh(2s) + cosh(2i))T (3) 
-4(sinh(2s) - sinh(2t))T (4) . 



(86) 



Note that the last form factor is antisymmetric and does not contribute at zero temper- 
ature. 

Now we consider II 12 and II 21 . From (Pill) we get 

n A 2 , + Ilfy = 5^(k - 2p) x G^(p)2 (5\'^k u i - <WV) G vv{p ~ k) 
+2 (5 Xtl K ~ huk^) G^{p)5^ v >{k - 2p)x>G v * v (p - k), 

and further 

K^+M 2 1,)9 T = (2{(k-2 P ) x ((E^-E- s - t )k) y 

+ - £r-*) T fc) (k - 2p(s)) v } . (87) 

We use the averages fl60l and obtain, 



+(^)4( i -vr- 2 " 



6r. 



This is an expression of the form of ( 1A.4I) with Pa from Eq.( !77l) and with additional 
M-dependent terms, where 

Q = E s+t - E- 3 - 1 = 2iF sinh(2g). (88) 
Then from Eq. (1A.6j) and accounting for the ^-dependent structures we find 



M 12 (s,t) + M 21 (s,t) = 4sinh(2g) 



~2aT^+rT^ 



+ (2(uk) - BJT^ -|t( q ) 



(89) 



Next we consider the contribution of II 33 together with the contribution from the 
ghosts, rF host . We get 



n 



33 



(<WG° ~~ + Sx'v'Pn') G u 'u(p - k). 



(90) 
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We use the property of the propagator 

P„G^(p) = G(p) Pll , (91) 
and obtain after simple calculation, using, for instance, the cyclic property of the trace, 

nfy = G X y(p)G(p-k)(p-k) 2 +p 2 G(p)Gy X (p-k) 

+P\G(p)(p - k) y G{p -k) + (p- k) x G{p) Py G{p - k). 

In the first two lines in the r.h.s. one line collapses into a point by means of, e.g., 
p 2 G(p) = 1 and the corresponding graph becomes a tadpole like contribution which will 
be considered below separately together with the tadpole diagram. The last two lines in 
the r.h.s. are just equal to the contribution from the ghosts, second line in (HUI) . with 
opposite sign and cancel. So we obtain 

II 33 + n ghost = 0. (92) 

Let us turn to the II 13 and II 31 contributions. We start from 

n Iv = - 2 P)\G^(p) (5av(p ~~ k )y' + <WP M ') G v i v {p - k), 

nfv = (5a m G° - k) v + S\ uPlx ) G^i{p)5^ V '{k - 2p)x'G u > u (p - k) (93) 

and arrive at 

MH(s,t) = -((k-2p) x (E»y(k-2p(s)) + E as k) x ,e) T , 
M&{*,t) = -((E s y(k-2p)-E as k) x (k-2p(s)) x ,Q) T , 

where the notation 

E sy = - (E s+t + E- s - 1 ) = <y" + 8 1 cosh(2g), 

= - (E s+t - E- s - 1 ) = iFsinh(2g) = -Q (94) 



E 



is introduced. Remind that the contributions appearing after the averaging procedure 
can be obtained by the substitution p^ — > p^ + = [(4fc) + u]^ in the final expressions. 
Doing so we obtain for the averages 



M A 13 ,(M) = l.-P x Q x , + Ai (^pZ ] ^0 



T- 



XX' 



M 3 x h(s,t) = {-q a p a T + \Q 



XX' 



T, 

(95) 



where Pa is given by Eq.( !77l) with the mentioned replacement being done, and Q x = 
Q x — 2u x with 

Q » =(( £M+£!y ( i -4) T )*) A , 

= rly + idy ^ sinh(2g) — a cosh(2g)^ + ft cosh(2q)hy 

= S ly + 7 idy +5 hy. (96) 
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The second contributions to the both lines in the r.h.s. in Eg. (1951) must be integrated by 
parts. 

This expression differs from that of in Eq. (f78|) by the factor E sy . That results in the fol- 
lowing replacements in the Eqs. fl7H|) . fl8TT) . (3 — > /3cosh(2g) = 5, — > c °^j j cosh(2g) = 
7. With these substitutions done and accounting for the overall factor 8, we calculate this 
contribution, 

Tfsy Zf—s rp 

8i — — = [2Bi - 4(uk)](r^ + 55 1 - - ijF). (97) 

Again, we have chosen the constant C = —2j^- and omitted the surface term. The latter 
will be considered separately. The above expression gives the parameters a,b,c entering 
the form (IA.4I) for this part. By using Eq. (jA.6j) from Appendix and the structure of the 
it-dependent tensors, we obtain 

n is + n 3i = 2r 2 T (i) 

+2 { \(3 2 - a 2 } cosh(2g) + a sinh(2g) } T (2) 

+r/3(l + cosh(2g)) T (3) 
-2r sinh(g) sinh(£)T {5) 
+ [2(uk)(a - 7 + 27) 



+/ 2 r( 7 - a - 27) + h 2 /3(>y - acosh(2g) - 27)] T (6) 
iN 
qf 

+r(a - 7 )T (4) + (7 - a)T (9) + fc 4 /3(l - cosh(2g))T (10) I . (9* 



2—T {7) + / 3 ( 1 + cosh ( 2 g)) T (8) 



Finally, we need II 23 and II 32 . Proceeding in the same way as before, we derive from (T4"0"j) 
and gl]) 

IT 23 , = 2{8xnK - 8 Xv k^)G^{p){8x'n'{p - k) v > + 5 X 'v'P^)G u i u (p - k), 

Uf x , = 2(5 XfJ ,(p-k) l , + 8xvP tl )G flfl >(p)(5x> fl >k v i -Sx'v'k^G^vip-k), (99) 

which gives 

M 2 x 3 As,t)(Q) = (2/ - E$ (kE\k-p{s))) - E{ y (kE- s p(s] 



T 



+ (E- S p(s)) x {E-'k)^ + (E t (k-p(s))) x (E s k) x , |6 
Ml 2 As,t)(e) = (2{ - E$ (kE-\k-p)) - E{ x , ( P E~ s k)) 

+ (E- s k) x (E- t (k-p)) xi + (E t k) x (E s p) xl \Q) . (100) 



Now the average (l60j) must be used. As it was noted, this results in the substitution 
Pn — ► P^ + Ufj, in Eq. fllOQI) . For the operator p(s) the matrix (^t/c) must be substituted. 
In this way the ^-dependent part is calculated. The w-independent part can be rearranged 
according to 

Mil + Mil = 2{A + B) (101) 
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with 



A 



B 



-2E{ X , (k^k) + (E*k) x UEt - 1 



D / A 

E s - 1 
-2E7Z ( k^—k 



+ 



E s -1 



E~ s k 



D 



E s - 1 



D 



k) (E s k) x , 



A and B have the structure of (IA.4I) . For A we obtain 

r' = 1, a' = sinh(2t), = cosh (2^ 

, t , 1 r/ cosh(s) sinh(t) 

s = ~, 7 = o«> 5 = ^7"^ 

q 2 sinn(g) 

and 



(102) 



a! = -2 [U 2 + 5'h 2 j , 

6' = cosh(2t)a', c = sinh(2g)a , 



(103) 



and for i? 



and 



1, a" = — sinh(2s), /3" = cosh(2s) 

s .. . .. sinh(s) cosh(t) 

-, 7 = -7,(5 = — — 

q sinn J 



(104) 



a" = -2\^l 2 + 5"h 2 j, 

b" = cosh(2s)a", c" = - sinh(2s)a". 



(105) 



For convenience, in the formulas f II 2 j) and fjl04j) we used apostrophes for denoting similar 
parameters having the same structure in different parts of Eq. (1A.4I) . 
Putting these contributions together we obtain with (jA.60 



n 23 + n 32 = 2 |_ 2T (l) 

cosh(2t) cosh(s) sinh(i) + cosh(2s)B cosh(t) 



sinh(g) 



+2 cosh(£)Bsinh(t) 



T (2) 



1 + 



s cosh(2s) + t cosh(2t) 



T (3) 
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s sinh(2s) — t sinh(2t) ^ 
Q 



+ 
+ 



■1 + 



s cosh(2s) + t cosh(2t) 



T (5) 



,ssinh(2s) — tsinh(2t), j9 n , . . 
— —)l 2 + 2 cosh(g) sinh(£)/i 



T ( 6 ) 



sm 



h(g)cosh(QT( 6 )- Sinh(g) ; Sinh(0 T( 8 ) 

K4 



+2(wfc) 

- sinh(g) cosh(£)T (9) + sinh(g) sinh(£)T (10) 



(106) 



Now we consider the contribution of the tadpole diagram, the last line in Eq. (l40l) . 
Accounting for the explicit form of the propagator, 



GM = J dqe-^E i 



q 



(107) 



and calculating trG(p) — 2 / dqe qp2 (1 + cosh(2g)) and the bracket averages, we obtain 



fjLV 



{G[iv + — 2trG5 ^J) T 

-A Art r r2 ' 



E 

N=—oo 



dqe ^ 



[6^(1 + 2cosh(2g)) + #,(2 + cosh(2g))] 
q sinh(g) 



(108) 



The contribution from the tadpole like terms coming from Ilj^, Eq. fl90l) equals just the 
first two terms in the Eq.( 1108p . That adds the extra terms 2(S^ U + 5^ u cosh(2g)) into the 
numerator of the above expression and cancels the first and the last terms in the total. 
Hence, the final contribution of the tadpoles and the tadpole like terms is 



1 +2? roo §\\ cos h(2g) + 5 J 



tot 



I f°° iV z 

4- — — V / dq e~^r 
(An) 2 „^„ Jo H 



N=-oo 



gsinh(g) 



(109) 



Now let us consider the surface contributions from the sum IP^ = IT 11 ) + IT 13 ) + n( 31 ); 



4 / dsdt 

surface 



^cosh(2 g )) + ^^ 
q sinn^gj 



cosh(£ ) u ■ u . . , 
-iF— fl-2-— cosh (2q) 
smh(g) (uk) 



(110) 



where the explicit expressions for the parameters r, (3 from Eq. (1771) are substituted. To 
relate the integrals in q and s, t we introduce an integration variable v = £/q in place of q. 



00 +1 

In this case J| = ^Jjj, and / dsdt — > | / gdg / dv. Hence, the integration over £ results 
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in an integration by parts over v. In the new variables the above equation reads, 



-%F 



dqQx 

cosh(gv) 
sinh(g) 



6hcosh(2q) + 5~ 



sinh(gw) 
sinh(g) 



2T^r-cosh(2g) 



{uk) 



with 



where 



e. 



exp 



1 



2(uk)-q(l-v) 



N 2 



(0) 



and 



ST 
S+T 



tanh[±q(l-v)] tanh[±q(l+v)] 
tanh[i<j(l-i>)]+tanh[ig(l+i))] ' 



AqT 2 

(47r) 2 gsinh(g) 
Hence it follows that 



(©), 



(112) 



(113) 



(6 



T 



1 



e iiT 2 



(4ir) 2 q sinh(g) 

| iN 



e 4< ? T ^ 



(An) 2 gsinh(g) 
After substitution of these functions in Eq. flllll) we obtain 



n 8 ' = 2 / dq (©T 

n N=-oo 



+ 1 



cosh(2g) + 5 J 



1 + e 



iNk^ 



(114) 



/ U • 11 

- iF coth(g) + 2—— cosh(2g) 



l-e 



(115) 



Now we note that in the imaginary time formalism the external momentum /c 4 = 27TOfcT 
with rik = 0,±1,±2, ... The phase factor in the exponential is = 2irNnk- So, the 
exponentials in the brackets equal to 1. This results in the factor 2 for the first two terms 
and zero for the third term in the integrand. To find the last term we also note that for 
k 4 ^ the denominator (uk) ^ and the bracket is zero too. For A; 4 = 0, we expand the 

Thus, we come to the 



exponential in a series and find — 2 1 — e t 



;A r A'4 



u-u 
(uk) I 



_ _N 2 
\k 4 =0 — -^fl 

conclusion that the last term contributes in the static limit k± — 0, only. 

As a final step, we take the Eq. (II 151) and the tadpole contribution Eq. (11091) together. 
The terms in front of 5 n and S x do cancel in the total. So, the only contribution is the 
M-dependent part coming from IF '* . We denote it as U s t 



Hot) 



Ly tot 



(4tt) 



u-u cosh(2o) sL 

d( l E T^tt ■ w e ^ 
N t^oo ( uk ) qsmh{q) 



iNk^ 

1 — e~^~ 



(116) 



This term contributes in the static limit only and it is transversal by itself. This is because 
being multiplying by k\ it is zero due to the difference in the curly brackets, if ^ 0. 
The integrand is nonzero if k A = 0, but for the product k\(Il^ t ) \y = holds. In fact, it is 
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the limiting expression of the tensor T^ 8 > for fc 4 = 0. However, we will consider this part 
separately for convenience. Thus, we collect all the contributions coming from individual 
parts. 

Now let us gather them together to obtain the one-loop form factors M^'(s, t) for the 
polarization tensor. We present the results as the list of explicit functions of variables 
q = s + t and £ = s — t: 



Mi 

M 2 

M 3 

M 4 
M 5 



4-2 



cosh(2g, 



1 — cosh(g) cosh(£) 
(sinh(g)) 2 

£sinh(£) 



-2 cosh(2g) 
2- ( sinh(2g) 



gsinh(g) 

cosh(g) — cosh(£) 



2 + 8 cosh(g) cosh(£), 
2 + 6 cosh(£) cosh(g), 

- 6 cosh(g) sinh(£), 



sinh(g) 



-2 + 2cosh(g) cosh(£), 



M, 



(i) 



coth(g)(l — 3(sinh(g)) 2 ) + sinh(£) cosh(g) 



+ 2 



Q 

sinh(£) 
sinh(g) 



I 1 



coth(g)(l — 3(sinh(g)) 2 ) + 2sinh(£) cosh(g) 



MP 
M 7 

M 9 



iN 

—k 4 2 (sinh(2g) - coth(g)) 



iN 1 C ( 
qT /c 4 q 



-2cosh(2g)), 



^1-2 
qT k 4 \ 

iN 



sinh(£) 



o. 



sinh(g) 
cosh(g) — cosh(£) 
sinh(g) 



4sinh(g) sinh(£)J , 

sinh(2g) — 2 sinh(g) cosh(£) 



(117) 



The symmetric form factors have to be multiplied by and the antisymmetric ones - 
by G^, when the integration over s, t is carried out. It is interesting that M\q is zero in 
one-loop order. 

Thus, according to Eq. fl38|) . we presented the polarization tensor in the form 

= ET«n«(A;) + (n i l); 



MA' 



1=1 



U (i \k) = dsdt M®(s,t)e T 



(118) 



N=-oo 



as double parametric integrals over the proper time parameters s, t and the temperature 
sum. The last term in the fist line is written in Eq. flll6j) . This representation is crucial 
for what follows. It preserves as much symmetries of the polarization tensor as possible. 

The obtained expression for the polarization tensor could be used in various applica- 
tions. In the next three sections we consider the zero field limit, the Debye mass and the 
magnetic mass of gluons in the magnetic background field at high temperature. 
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6 Limit of zero background field 



To make a link of our formalism with a usual one, let us consider the limit of zero 
background field, B = 0. In our dimension less variables this simply corresponds to the 
limit q, £ go to zero. In this case, the most form factors and operators also go to zero. 
More precise, we have for form factors, 



.I/, M 2 M 3 1 



2 r 



(119) 



and 



M 7 = M 8 



iN 1 
qT fc 4 



9, 



(120) 



All the other operators or form factors are zero. In accordance with Eqs. (1351) . 0371) this 
means that we obtain two transversal operators with the correct formfactors at finite 
temperature [15]. Hence, at B = 0, the polarization tensor can be written as 




dsdt 



e i (2{uk)t — Kf) 



iN 1 ^ 



K xy -2[^^\B 



qT fc 4 q 



(47T) 



N=—oa n v 



u\U\i e 4< ? T 



(uk) q 2 



XX' 



;m) 



where we substituted the surface contribution Eq. (jll6j) at B = 0. 

In the representation in terms of a series resummed according to Ref. 
N = corresponds to the zero temperature case. In Eq. (11211) this is 



the value 



n 



(T=0) 
AA' 



dsdt 



_ k 2st 

e i 



4 




(122) 



that coincides with Eq.(103) in Ref. [10]. This part must be renormalized in a standard 
way. Actually, normalizing to B = at T = it must be simplify skipped. 

Now, let us calculate the Debye mass squared defined as the limit of the form, 
m 2 D = — n 44 (T, fc 4 = 0, k — > 0)[15j. Within the representation (11211) . only the last term 
contributes and we obtain 



m 



D 



1 +2? N 2 
TZ2 1^ f 2 

^ n N=l 1 



dq 



e 4 " t2 



(123) 



The integral is simply calculated, / dq 

o 

J2 



AT 

4qT 2 



16T 4 



^3 — — -jfr, and the sum is expressed through 

Riemann's Zeta-function, £(2) = ^r. Thus, we obtain m 2 D = |T 2 , which is the well known 
result [15]. 

The next important parameter is the "magnetic" mass squared which can be deter- 
mined as the limit for transversal with respect to the external field direction components, 
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m magn = — H 12 (T, fc 4 = 0, k — » 0). In this case, the operator Ki 2 in Eg. (11211) contributes 
(the component Byi = 0). To calculate the form factor in the high temperature limit, it 
is convenient to make an inverse resummation according to Eq. (l57| . In the static limit, 
&; 4 = 0, the parameter a = and we have 

+°° w 2 s / n \ 1/2 +00 

^ e (-^) =27rT ^y £ e <r4**M*) m (124) 



AT=-oo v7F/ N=-oo 



For the components II 12 we then get 



n.^)^Jj^/<M2,T)(iV /2 





2^ 



e 



(-*V(i-«)) e (-^4^T^) ( 4 - 2 m I tf 12 , (125) 



where new variables, s, t — ► s = git, t = q(l — u), were introduced. The high temperature 
limit corresponds to the value of N = in the Matsubara sum. In that case the integrals 
can be calculated easily. First we compute the integral over q which delivers T-function, 
r(|)/(fc 2 w(l — u)) l l 2 . Then the integration over u gives 3ir. So, for the form factor we 

obtain U(k)i2 = ff-- Hence for the Green function, 

follows. Here k is the length of the three-momentum vector k. This expression has a 
fictitious pole that was an old problem of gauge theories at finite temperature [15J. 

From the above calculations we see how to proceed in the present formalism. The 
procedure remains actually the same when the field is switched on. Bellow, we investigate 
an influence of the field on the Debye mass and the existence of the fictitious pole. 



7 Debye mass in the presence of the background field 

The gluon Debye mass squared in the background field is defined as before, m 2 D = 
— n 44 (T, B, fc 4 = 0, k — > 0), where we have to use the expression Eq. fl!16j) from the 
polarization tensor. In this limit it reads 

= 55/ j E ^ — ^f* (127) 

where the dimensional parameters are restored. We investigate the case for the field and 
temperature ratio s = <C 1 and present Eq. (1127|) in the form: m 2 D = |T 2 /(s), where 
the function f(s) is to be computed. It satisfies the condition /(0) = 1. With this 
parameter introduced we return again to the dimension less variable Bq — > q and write 
f(s) in the form, 



' 1 ^ r l cosh(2g) 1 



71 J Q H N=l 



q sinh(g) q 



(128) 



25 



where the first term in the square bracket delivers the zero field limit /(0) = 1. To 
calculate the other terms resummations will be done. By differentiating Eq. (l57I) with 
respect to Z with a = 0, Z = -, we rewrite f(s) as follows 



1 3 1 

1 + 

Z 7T Z 



+ 



H + f k - 

/ iV=-oo 



sinh(g) q 



1 



7r 2 ;vy 



(129) 



Here we made a rearrangement of the integrand in the above equation and split it into 
three parts - the tachyonic one ft{s) coming from e q , f'2{s), and fs(s) coming from the 
curly brackets. Then we calculate the leading terms of the high temperature expansion 
which is given by the terms with N = in the temperature sum, next to leading - by the 
terms with N > 1. 

Now, let us calculate contributions from N = 0. For f t (s, N = 0) we get 



Z 7T Z 



dq 



(130) 



This integral diverges at the upper limit that reflects the tachyonic mode in the spectrum 
of charged gluons (I2l)j) . To obtain its value we make the inverse Wick rotation in the 
g-plain, that is, replace q — > qe~ l7T . After that a simple integration yields 



The second term is 



ft(s,N = 0) 



f2(s,N = 0) 



3 1 3 
-i-^y/Ws^ = -i— Vs. 
2 7i z Zir 



3 1 
~2V 2 



and can be easily computed, f2(s, N = 0) 



2tt 



dq 



V'i. The third term is 



(131) 



(132) 



/a (s,N = =--v™/ — - 
2 7H J , A? \ smh 



y/q ysinh(g) g / 



that can be expressed trough a Zeta- function, / 3 (s, iV = 0) = |^ 
we derive 



(133) 



7rsV27rC(i|)- Finally 



fs{s,N = 0) 



V2n 



Similar simple integrations for the value N — 1 result in 

n 25C(3) 



Thus, for the function f(s) we obtain 
f(s)= 1 + 
+ 



4 ?H 



s 2 . 



(134) 



(135) 



V2 



71 



3 

2t1 



25C(3) 9 3 
s — z- 



4 7T 4 



2tt 



(136) 
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Hence, for the Debye mass we derive 



m 2 D {B) 



3 



1 - 0.8859 



2T 



0.4775 



B 2 



V16T 4 



-i 0.4775 



(#) 


B 3 1 







(137) 



where the numeric values of the coefficients in Eg. (11361) are substituted. This expression is 
interesting in two respects. In presence of the field, the screening mass for color Coulomb 
forces is decreased as compared to the zero field case. This behavior was already deter- 
mined in Ref. [7j in calculation of other type. The imaginary part of m 2 D is signaling the 
instability of the states because of the tachyonic mode in the spectrum (1251) . However, 
the numeric value of it is small as compared to the real one. It is of order of the usual 
plasmon damping factor at finite temperature. 



8 Transversal modes in the field presence 

Let us investigate the behaviour of the transversal modes in the field at high temperature. 
We have to calculate the mean value of the polarization tensor in the states given by 
Eq. (l26l) for the polarizations s = 1 and s = 2 in the limit of = 0, k — > 0. Accounting 
for Eqs. (l3"41) . (l3"oT) . (I2"6"j) . we derive for the mean values 

(s = i\n{k)\s = i) = h 2 u 2 + z 2 (n 3 - n 5 ), (138) 



<* = 2|n(*)| a = 2) = ^in x + (h 2 + ^A 2 ) n 3 

+ —gr*-^ ( 139 ) 

where U = k±, k 2 = h 2 + I 2 , 1 2 = l 2 + l\. To consider the behaviour of the static modes in 
the perpendicular with respect to the field plane we put h = and A; 4 = and get 

(s = l|n(Jfe)|s = l) = h 2 Tl 2} 

(s = 2|n(A;)|s = 2> = h 2 (U 3 + U 5 ). (140) 

We have to calculate the form factors n 2 ,n 3 and n 5 educed in Eqs. (lll7j) . (lll8p also for 
this case. 

The procedure of calculations is quite similar to that of in the previous sections. We 
describe its steps by computing the form factor n 5 . To carry out integration over £, we 
change variables, s = qu, t = q(l —u), as in section 6. In the limit of interest, h 2 —>■ 0, we 
restrict ourselves to the leading in this parameter term when the function (0) Eq. (1461) is 
substituted. More precise, we take into account the first term in the expansion, 

(0) 1 i + OM (Ul) 
(47i) z gsmh(g) 

because, being substituted into Eq.( ll40l) . the ~ 0{h 2 ) terms result in a next to leading 
correction. Then we make resummation of the series in iV and take into consideration the 
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term with N — 0, that gives the leading high temperature contribution. In this limit we 
obtain for the form factors 



1 oo 

U {N=0) = -J — du I q ^ MA a u) (142) 

(47r) 3 / 2 Vs/ J sinh(g) Ug ' Mj ' 1 j 

where the functions from Eq. (II 171) should be substituted. Here we restored the overall 
factor -j=. All other variables are dimension less. 

For the function M 5 , integration over u results in the expression 

M 5 (q) = -2 + 2 cosh (?)sinh(g) (143) 

1 

which after substitution in Eq. ( 1142ft leads to the integral 

dq 



V + — + 



sinh(g) Jq Jq_ 



(144) 



The second and the third terms are calculated in Eqs . (1 1 30 1) - fl 1 32 p in section 7, 1^ = y/n, 



I5 = —i\^TT, and for the first one we compute 1^ = ~(— 4 + \/2) Jlx C(|)- Thus, we 
obtain in the total, 

IT 5 = j^y- 2 ^ [-4.21405 - 1.77245.] , (145) 

where numeric values of the integrals are substituted. 

Similar calculations carried out for the form factors H2 and H3 yield 

II 2 = ^-^-^[-5.79894- 7.08982.], 

n 3 = ^L_-^ [1.04427 -8.86227,]. (146) 



The above expressions have to be used in Eq. (I140j) to obtain final result. The sum of 
n 3 + n 5 equals, II3 + II5 = ^= [—3.16978 — 10.6347z]. The imaginary part is signaling the 
instability of the state because of the tachyonic mode presenting in the vacuum, and the 
real part is responsible for the screening of transversal gluon fields. 

As concerns the imaginary part. Since its origin is clear, one should conclude that 
the one-loop result needs in additional terms coming from some resummation at finite 
temperature in the sector of charged components similar to that of in Refs.[7],[8]. That 
we left for the future. 

Let us turn to the real part and substitute it in the Schwinger-Dyson equation 

D-\k 2 ) = k 2 -U(k) (147) 
for the neural gluon Green function. We obtain for the mean values 
( s = 1 \D-\h 2 )\ s = 1 ) = h 2 -Re(Tl 2 )h 2 

= h 2 (l + 5.79894-^=1 (148) 
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and 



( S = 2 \D-\h 2 )\ s = 2) 



h 2 - Re(U 3 + n 5 ) h 2 
h 2 ( 1 + 3.16978-?=] 




(149) 



These are the expressions of interest. 

Two important conclusions follow from Eqs. (ll48j) . (11490 . First, for the transversal 
modes in the field presence, there is no fictitious pole similar to that of in Eg. (11261) . 
Second, there is no the magnetic screening mass in one-loop order. The transversal 
components of gluon field remain long range in this approximation, as at zero external 
field [15J. These conclusions are in a complete correspondence with the results obtained in 
Ref . [14j in the Fujikawa gauge. In this paper only the static limit = was investigated. 
As it has been discovered in Refs.[7J,[B], for the charged components the magnetic mass in 
the color magnetic field at high temperature is generated in one-loop order. So, there is 
an essential difference in this point. It is also interesting that the non transversal tensor 
structure contributes to the transversal state s = 2. 

9 Discussion 

In the present paper the operator structure of the neutral gluon polarization tensor in 
a constant Abelian chromomagnetic field at finite temperature is derived. There are 
ten operator structures contributing in general case. We mention that these structures 
do not depend on the gauge group; for a SU(3) these are the same as for the SU(2) 
considered here. As investigated in details in Refs. [I0],[TT] at zero temperature, the 
polarization tensor is not transversal. The tensors T^\T^\T^ and satisfy the 

weak transversality condition Eq. (l32p . As actual one-loop calculations showed, the form 
factors 115,116,119 are nonzero and II10 = 0. So, only three non transversal operators 
contribute in the total in this approximation. 

As at zero temperature Ref.jlO], to carry out integration over internal momentum 
of diagrams in the field presence which is necessary in calculation of form factors, the 
Schwinger algebraic method [T2] was applied. It has been extended to account for finite 
temperature. As a result, we have derived the explicit expressions for integrals which 
present the form factors in terms of two parametric integrals over the proper-time param- 
eters s, t, as at zero temperature, and extra series in discrete variable N accounting for 
the temperature dependence. This representation preserves gauge invariance and even 
Lorentz covariance at each step of calculations as far as possible. It is simple and conve- 
nient for investigations. 

Within this representation, the one-loop form factors for the polarization tensor of 
neutral gluons have been obtained in the background Lorentz-Feynman gauge. Note here 
that in previous papers [17] . [IB], [2] this tensor in the field without and with temperature 
included has been computed in the Fujikawa gauge. In this non-linear gauge the structure 
of the neutral polarization tensor is simpler as compared to the Lorentz gauge and the 
tensor is transversal, but the charged sector has a much more complicated structure. 

The tensors Eqs.( l3"4"]) . (l36"I) and form factors (I117p . (lll8p present the structure of the 
neutral polarization tensor off shell and can be used in various applications. First of all, 
the spectra of gluon modes in the field at high temperature can be investigated. They 



29 



may also serve as the basis for different kinds of resummations of perturbation series. 
This is because to make a resummation one needs in the Green functions of shell. As 
far as we know, the operator structure of the gluon tensor in the considered background 
was not educed before. To relate our approach with a common one at finite temperature, 
we have investigated in short the zero field case and recalculated the corresponding well 
known results. 

As applications, in the present paper we restricted ourselves to the consideration of 
two important parameters - the Debye mass of neutral gluons, Eq. (11371) . and the magnetic 
mass of transversal neutral gluons, Eqs. (11481) . (I149j) . in the field at high temperature, in 
one-loop approximation. It was found that all the asymptotic expansions for the high 
temperature limit are expressed in terms of Riemann's (^-function. As it follows from 
Eq. fll37p . the screening temperature mass of plasmons is decreased as compared to the 
zero field case. That increases a radius of color Coulomb forces when the background 
magnetic field is present. This behavior may be important for quark-gluon plasma. As 
concerns the magnetic mass, from Eqs. (I148p . (11491) it follows that it is zero in one- loop 
order, as at zero field. At the same time, the fictitious pole, appearing in one-loop order 
at B = (see Eq. (ll26l) ). disappears in the field presence. This is important property of 
hot gluon plasma. Already earlier one believed that the elimination of the fictitious pole 
requires a resummations of perturbation series [IS], [IE] and it is removed by the gluon 
magnetic mass. 

In connection with the last result we would like to speculate a little. If one takes the 
one-loop result seriously and the magnetic mass is zero, that means that neutral gluon 
fields remain long range like Abelian magnetic fields at high temperature. This is in 
contrast to the charged gluon fields which acquire the magnetic mass in the field at high 
temperature in one-loop order [Zj,[ll] and therefore are screened at long distances. This 
behaviour is important either for problems of quark gluon plasma or the early universe 
because at high temperature the stable Abelian magnetic fields are spontaneously created, 
as it has been derived in the daisy resummations in Refs. [5] , [6] , [7] , [8] . 

To make a final conclusion about this phenomenon, one has to apply the super daisy 
resummations which can be realized on the base of solution of the Schwinger-Dyson equa- 
tions for the common system of the neutral and charged gluon Green's functions. The 
general structure of the Green function for neutral gluon fields is derived in the present 
paper. The case of the charged gluon field will be considered separately elsewhere. 
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Appendix 

In this appendix we collect formulas which are used to identify the contributions to the 
form factors, i.e., the contributions which go with the tensor structures TffJ. 

In the course of calculation, from the graph (Fig. [2]) contributions appear which have 
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up to a constant the following structure. First, 

T{ x , = P\P\i + a4, + b8& + icF xx , 7 (A.l) 

where Pa is given in terms of the vectors (|33|) . 

P A = r l x + a id x + (3 h x (A.2) 

and r, a and (3 are some functions of the variables s and t. The transposition in P x , 
changes the sign of d\>, P x , = r l x — a id\ + (3 h x . 

It can be seen that the expression in Eq.QSBH fulfills ([32]) if (rl 2 + f3h 2 f + a/ 2 + bh 2 = 
holds. In that case it can be represented in terms of form factors according to 

Fly = -r 2 T« + (a 2 - i3 2 )TS» - r^f, - raTg! 
+ r(rP + jW) + a r ,g + + + (g 

which can be checked by inserting the explicit expressions (11031) . 

A second type of expressions appears which has a slightly more complicated form, 

•P\\> = PxQx' + QjPy + aSlx + b5i x , + icF xx , (A.4) 
with P\ from (IA.20 and 

Qx = s l x + l idx + 5 h x , (A.5) 

and s, 7 and 6 are also some functions of the variables s and t. In parallel to the above case, 
if for (JXIU) the condition ([32]) is fulfilled than (a+2rsl 2 )l 2 + (b+2(35h 2 )h 2 + (r5+s(3)2l 2 h 2 = 



must hold. In that case the representation in terms of form factors is 

T 2 xxt = -2rsT$, - 2((35 + ai )T$, - (rS + sP)T$, 

+ (r 7 - sa)T$ + ((a + 2rsl 2 )^ + r6 + s/?) rg 
+ (c - (r 7 - sa)/ 2 + (a<5 - /3 7 )/i 2 ) P? A ,. (A.6) 
Formulas (1A.3j) and (1A.6[) are used in the section 5 for calculation of form factors. 
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